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Abstract 

The BAEL4R detector at the PEP-II asymmetric B Factory at SLAC collected a sample of 
23 ■ 10 6 BB pairs in the years 1999 and 2000. Using this data sample, we measure the amplitude 
of the time-dependent CP-violating asymmetry in neutral B decays to the CP eigenstates J/ipK®, 
^{2S)K° S and J/ipK^. We find a value of sin2/3 = 0.34 ± 0.20(stat) ± Om(syst). We also present 
preliminary measurements of the B°B° oscillation frequency and of the lifetimes of charged and 
neutral B mesons. 
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The measurement of CP-violating asymmetries in the time distribution of decays of neutral 
B mesons can provide a direct test of the standard model of electroweak interactions § and is the 
primary goal of the BABAR experiment at the PEP-II asymmetric-energy e + e~ collider at SLAC. 
Decays of B° and B° mesons into charmonium CP eigenstates due to b — * ccs transitions (e.g. 
B° — > J/ip Kg) can be used to measure sin2/3 (where (3 is one of the interior angles of the unitarity 
triangle) with negligible corrections fromstrong interactions. 

In the years 1999 and 2000, BABAR ! collected a sample of 23 • 10 6 BB pairs (20.7 ftT 1 on- 
peak) at PEP-II, where BB mesons are produced at the T(4S) resonance in collisions of 9.0 GeV 
electrons and 3.1 GeV positrons. The boost (f3j) = 0.56 along the collision axis (z) resulting from 
the asymmetric beam energies allows the determination of the proper decay time difference At of 
the two B mesons from the measurement of the decay length difference Az, whose average value is 
(fij) ct b o ~ 260 fj,m. 

In T(AS) decays, B°B° pairs are produced in a P-wave state and evolve coherently until one of 
the B meson decays. At that time (At = 0) the other B meson has the opposite flavor. In events 
where one of the B mesons, Bqp, decays into a charmonium CP eigenstate and the other, Btag-, 
decays such that its flavor can be determined, the expected decay-time distribution (J 7 -) for 
events where the flavor tag is a B° (B°) is given by 

JF ± (At rec ;r, Am d ,w,sin2(3) = ^Ve' 1 *^ [1 TVCp(l ~ 2u>) sin 2/3 sin Am d At] ®K(6 t ;a) . (1) 

The mistag rate w is the probability to wrongly determine the flavor of Btag- The decay-time 
distribution is convoluted with a time resolution function TZ(St = At rec — At; a) with parameters a 
in order to account for the finite resolution of the detector, rjcp is the CP eigenvalue of the final 
state (-1 for decay modes with K®, +1 for modes with K®). 

For events where one of the B mesons decays into a fully reconstructed flavor eigenstate Bg av , 
the decay-time distribution for unmixed (TC+) and mixed (H-) signal events is 

H±(At vec ;T,Am d ,w,a) = -re~ r|At| [1 ± (1 - 2w) cos Am d At] ® K(8 t ; a) , (2) 

where an unmixed event is one where the Bfl av and Bt ag mesons have opposite flavor. The mistag 
rate and resolution function (the latter is dominated by the reconstruction of the Bt ag vertex) are 
the same as for the Bqp sample, and both can be determined from the Bfj av sample. 

CP candidates are reconstructed in the decay modes J/tftK®, Tp(2S)K®, and J/tpK® and are 
required to have a difference AE between the energy of the Bqp candidate and the beam energy 
in the center-of-mass frame of less than 3 standard deviations from zero. In addition, modes 
with a K® must have a beam-energy substituted mass tties = \J (-E'boam) 2 — {p c pt) 2 > 5.2 GeV/ c 2 
(tkes > 5.27 GeV/ c 2 for candidates counted as signal). The distributions of tues and AE for the 
Bqp candidates are shown in figure p]. A sample of B decays, Bfj av , reconstructed in the flavor 
eigenstate modes f\ h + (h + = ir + , p + ,af) and J/ipK*° (K*° — > K + ir~) is used to measure the 
B° lifetime and Am d . A sample of charged B decays, B c h, in the final states J/ip , il>(2S)K + 
and D^°tt + is used to measure the B + lifetime and serves as a control sample. Yields and purities 
for events with a flavor tag are summarized in table 1. 

The vertex of the other B in the event (Btag) is determined by fitting the tracks not belonging 
to the reconstructed Bqp, Bg av or B c h to a common vertex. Tracks from 7 conversion are removed 
and reconstructed K® and A candidates are used as input to the fit in place of their daughters. 

a Throughout this paper, flavor-eigenstate decay modes imply also their charge conjugate. 
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Figure 1: m^s and A_B distribution for Bqp 
candidates with (top) and Kl (bottom). 



Sample 


^tag 


Purity (%) 


sin2/3 


Full CP sample 


529 


69 ±2 


0.34 ± 0.20 




256 


39 ±6 


0.87 ± 0.51 


Jft>K», *P(2S)K" S 


273 


96 ± 1 


0.25 ± 0.22 


- Lepton tags 


34 


99 ±2 


0.07 ± 0.43 


- Kaon tags 


156 


96 ±2 


0.40 ± 0.29 


- NT1 tags 


28 


97 ±3 


-0.03 ± 0.67 


- NT2 tags 


55 


96 ±3 


0.09 ± 0.76 


- B° tags 


141 


96 ±2 


0.24 ± 0.31 


- B a tags 


132 


97 ±2 


0.25 ± 0.30 


-Bflav sample 


4637 


86 ±1 


0.03 ± 0.05 


Charged B sample 


5165 


90 ±1 


0.02 ± 0.05 



Table 1: Number of tagged events, signal purity, and result of 
fitting for CP asymmetries in the full Bcp sample, in various 
subsamples, and in the Bfl av and charged B samples. 



Tracks with a large (> 6) x 2 contribution are removed to reduce the bias from charm decays. We 
require a(Az) < 400 /mi and \Az\ < 3mm. The average resolution of Az is 190 /mi. 

To a very good approximation At ~ Az/c(f3j), but event-by-event corrections are made 
for the direction of the B with respect to the z direction in the T(4S) frame. The resolution 
function 1Z(5t;a) is parameterized either as the sum of three Gaussian distributions (sin2/3 and 
Awid measurements) or as the sum of a zero-mean Gaussian distribution and its convolution with a 
decaying exponential (lifetime measurements). In both cases, the event-by-event errors calculated 
by the vertex fits are used to scale some of the contributions to lZ(5 t ;a). From an unbinned 
maximum likelihood fit to the At rec distribution in the -Bfl av and B c ^ samples, including also events 
without a flavor tag, we obtain the preliminary results 

t b o = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps (3) 
t b + = 1.673 ± 0.032 (stat) ± 0.022 (syst) ps (4) 
t b+ /t b o = 1.082 ± 0.026 (stat) ± 0.011 (syst) . (5) 

The fit takes into account contributions from signal, background and outliers. The probability of 
each candidate to be signal is determined from the toes distribution. An empirical description 
based on toes sidebands and including both prompt and lifetime components is used to describe 
the At shape of the combinatoric background from other B decays and from continuum events. 

In order to determine the Stag's flavor tag for the sin2/3 and Arrid measurements, we use the 
flavor information carried by the charge of high momentum leptons (e, /i) from semileptonic B 
decays, of kaons from b — > c — > s transitions, of soft pions from D* decays and of high momentum 
charged particles not coming from the reconstructed Bcp or Sfl av candidate. Each event is assigned 
to one of four hierarchical mutually exclusive tagging categories (or else not assigned a flavor tag) . 
The Lepton and Kaon categories are characterized by the presence of an electron or muon with 
a center-of-mass momentum p* e > l.OGeV/c or p* > LlGeV/c, and of one or more kaons with 
a non-zero charge sum, respectively. The remaining two categories, NT1 and NT2, are based on 
the output of a neural network algorithm whose performance relies primarily on soft pions and 
on recovering isolated electrons and muons from semileptonic B decays. The tagging performance 
measured on data is shown in table 2. 

Based on the flavor of the reconstructed -Bfl av candidate and the flavor tag, events in the -Bfj av 



Category 


e (%) 


w (%) 


Aw (%) 


Q (%) 


Lepton 


10.9 ±0.4 


11.6 ±2.0 


3.1 ±3.1 


6.4 ±0.7 


Kaon 


36.5 ±0.7 


17.1 ±1.3 


-1.9 ± 1.9 


15.8 ± 1.3 


NT1 


7.7 ±0.4 


21.2 ±2.9 


7.8 ±4.2 


2.6 ±0.5 


NT2 


13.7 ±0.5 


31.7±2.6 


-4.7 ±3.5 


1.8 ±0.5 


All 


68.9 ± 1.0 






26.7 ± 1.6 



Table 2: Average mistag fractions Wi and mistag differences Awi = 
Wi(B°) — Wi(B°) extracted for each tagging category i from the 
maximum- likelihood fit to the time distribution in the Bfl av +Bcp 
sample. The figure of merit for tagging is Qi — £i(l — 2wi) 2 , where 
Si is the fraction of events in the i th category. The statistical error 
on sin2/3 is proportional to where Q = ^ Qi- 




At (ps) 

Figure 2: Raw asymmetry in the number 
of B° and B° tags for J/ip K° s , ip(2S) K° s 
(top) and J/ip Kl modes (bottom). The 
solid curve is the sin2/3 fit to these samples. 



sample are classified as mixed or unmixed. The time-dependent mixing asymmetry ^4(/\£ rec ) — 
{Nunmixed- ^ mixed) I (N unm ixed+ N mixed ) is shown in figure H Am d is determined from an unbinned 
maximum likelihood fit in which the mistag fractions Wi and Awi (8 parameters), signal resolution 
function parameters (9 parameters) and the fractions, lifetimes, dilutions and resolution function 
parameters of different background components (16 parameters) are determined simultaneously 
with Am d . The correlation between Am d and Wj is small because the latter are determined by 
events at low values of At where the mixing probability is small. 

An alternative method for measuring Am d is to use inclusively reconstructed dilepton events, 
i.e. events where both B mesons decay semileptonically and the flavor of each B is given by the 
charge of the high momentum electron or muon produced in its decay. Because of the relatively large 
semileptonic branching ratio and the high lepton identification efficiency, this method is statistically 
more powerful. The non-negligible backgrounds due to leptons from charm decays are minimized 
with a neural network technique which uses the lepton momenta and opening angle, and the total 
and missing energy as input. The resulting dilepton sample has about equal contributions from 
neutral and charged B mesons, but the former can be enhanced by an inclusive reconstruction of 
B° — > D* + i~v decays. The mixing asymmetry obtained with this technique is shown in figure ||. 

The preliminary results obtained for Am d with the two methods are: 

Am d = 0.519 ± 0.020 (stat) ± 0.016 (syst) % ps" 1 (B flav sample) (6) 

Am d = 0.499 ± 0.010 (stat) ± 0.012 (syst) % ps" 1 (dilepton sample) (7) 

The sin2/3 measurement is made with an unbinned maximum likelihood fit to the At rec distri- 
bution of the tagged candidates from the combined Bqp and B^ av samples with parameters similar 
to the ones used in the fit for Am^. The values of the B° lifetime and Amd are fixed to their world 
average a We find a value of □ 

sin2/3 = 0.34 ± 0.20 (stat) ± 0.05 (syst) . (8) 

The raw asymmetry in the number of B° and B° tags as a function of At rec is shown in figure |2| 
and has, as expected from eq. |], the opposite sign for CP even and CP odd modes. 

The determination of the mistag rates and At resolution function function is dominated by the 
high statistics -Bfl av sample. The largest correlation between sin2/3 and any linear combination of 
the other 34 free parameters is 0.076. The dominant sources of systematic error are the assumed 
parameterization of the At resolution function, due in part to residual uncertainties in the alignment 
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Figure 3: Mixing asymmetry A(At lec ) in the dilepton (left) and the Ba av (right) sample (preliminary). 



of the silicon vertex tracker, and uncertainties in the level, composition and CP asymmetry of the 
background in the selected CP events. The large Bcp sample allows a number of consistency 
checks, including fits to subsamples of the data by decay mode, tagging category and flavor of the 
Btag as shown in table 1. No statistically significant asymmetry is found in fits to control samples 
where no asymmetry is expected. 

The measured value of sin2/3 is consistent with the range implied by measurements and the- 
oretical estimates of the magnitudes of CKM matrix elements!. It is also consistent with no CP 
asymmetry at the 1.7a level. 
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